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ABSTRACT 

We present a theoretical study on the metallicity dependence of the initial—to—final 
mass relation and its influence on white dwarf age determinations. We compute a 
grid of evolutionary sequences from the main sequence to ^ 3 000 K on the white 
dwarf cooling curve, passing through all intermediate stages. During the thermally- 
pulsing asymptotic giant branch no third dredge-up episodes are considered and thus 
the photospheric C/0 ratio is below unity for sequences with metallicities larger than 
Z — 0.0001. We consider initial metallicities from Z = 0.0001 to Z — 0.04, accounting 
for stellar populations in the galactic disk and halo, with initial masses below ^ 3Mq. 
We found a clear dependence of the shape of the initial—to—final mass relation with 
the progenitor metallicity, where metal rich progenitors result in less massive white 
dwarf remnants, due to an enhancement of the mass loss rates associated to high 
metallicity values. By comparing our theoretical computations with semi empirical 
data from globular and old open clusters, we found that the observed intrinsic mass 
spread can be accounted for by a set of initial—to—final mass relations characterized 
by different metallicity values. Also, we confirm that the lifetime spent before the 
white dwarf stage increases with metallicity. Finally, we estimate the mean mass at 
the top of the white dwarf cooling curve for three globular clusters NGC 6397, M4 
and 47 Tuc, around 0.53Mq, characteristic of old stellar populations. However, we 
found different values for the progenitor mass, lower for the metal poor cluster, NGC 
6397, and larger for the younger and metal rich cluster 47 Tuc, as expected from the 
metallicity dependence of the initial—to—final mass relation. 
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1 INTRODUCTION 

White dwarf stars are the most common stellar evolution¬ 
ary endpoint. All stars with initial masses up to ~ IO.SMq 
(Smartt 2009, Doherty et al. 2014), more than 97% of the 
stars in the Galaxy, end their lives as white dwarf stars. 
Their evolution can be described as a simple and slow cool¬ 
ing process. Because white dwarf star are abundant and 
long-lived objects they convey important information about 
the properties of all galactic populations (Isern et al. 2001; 
Liebertet al. 2005; Bono et al. 2013). In particular, white 
dwarf stars can be considered as reliable cosmic clocks to 
infer the age of a wide variety of stellar populations, such 
as the Galactic disks and halo (Winget et al. 1987; Garcfa- 
Berro et al. 1988ab, Isern et al 1998; Torres et al. 2002), and 
globular and open clusters (Kalirai et al. 2001; Hansen et al. 
2002, 2007; Kalirai 2013). 

* E-mail: alejandra.romero@ufrgs.br 


Methods for determining stellar population ages from 
their white dwarf cooling sequences are usually based on the 
comparison of either the observed white dwarf luminosity 
function (Winget et al. 1987, Bedin et al. 2010, Garci'a-Berro 
et al. 2014), or the distribution in the Color-Magnitude Di¬ 
agram, with theoretical computations (see e.g. Hansen et al. 
2007). Both techniques rely on the use of extensive grids of 
white dwarf cooling sequences. However, the evolutionary 
stages prior to the final cooling curve must be taken into ac¬ 
count. In particular, a connection between the properties of 
the white dwarfs with those of their main sequence progen¬ 
itors is the initial-to-final mass relation (IFMR). This rela¬ 
tion can be applied to study the chemical evolution of galax¬ 
ies, including the enrichment of the interstellar medium, and 
the chemical evolution history of stellar populations in gen¬ 
eral. At the high mass end, the white dwarf IFMR can lead 
to constraints on the critical mass that separates white dwarf 
production from Type H supernova explosions, and at the 
low mass end, the relation represents a tool to probe the pro- 
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genitor properties of most of the evolved stars in old stellar 
populations, the majority of which are now low-mass white 
dwarf stars. 

The study of the IFMR started with the work of Wei- 
demann (1977) who compared theoretical models of mass 
loss to the observed masses of a few white dwarfs in nearby 
stellar clusters. This semi-empirical approach was followed 
by other authors. More recently, Catalan et al. (2008a) ob¬ 
tained a semi-empirical IFMR by re-evaluating the available 
data, that includes observations from individual stars from 
open clusters and common proper motion pairs. Kalirai et 
al. (2008) increased the observations and included data from 
two open clusters older than 1 Gyr and from the old open 
cluster NGC 6791, extending the initial mass range down to 
1.16M0. Finally, Kalirai (2013) presented an IFMR based 
on spectroscopic observations of white dwarfs belonging to 
open and globular clusters and Sirius B (see e.g. Kalirai et 
al. 2007, 2008, 2009), with initial mass ranging from 0.8M© 
to ~ GMq. In particular, Kalirai et al. (2008) bin the re¬ 
lation so that each star cluster was represented by a single 
point. 

In the last few years, the amount of data available has 
increased considerably. In the Galactic held, the Sloan Dig¬ 
ital Sky Survey increased the total number of spectroscopi¬ 
cally conhrmed white dwarfs stars to ~ 30 000 (Kleinmann 
et al. 2013, Kepler et al. 2014). In the Galactic Halo, recent 
Hubble Space Telescope observations of the globular clus¬ 
ters M4 (Richer et al. 2004; Hansen et al. 2004; Bedin et al. 
2009), uj Gen (Monelli et al. 2005), NGG 6397 (Richer et al. 
2006; Hansen et al. 2007) and 47 Tuc (Hansen et al. 2013) 
have similarly uncovered several hundreds of cluster white 
dwarfs. In all cases, the white dwarf samples are dominated 
by low mass stars. With the increase of the available data, 
the general trend in the empirical IFMR remains. However, 
the scatter in the data has increased, possibly due to a cor¬ 
relation between the mass loss rates and the properties of 
the host environment. Although this scatter could be re¬ 
lated to other properties as rotation, binary evolution and 
the presence of magnetic helds (Weidemann 2000), it is more 
likely produced by differences in the chemical composition, 
i.e. different metallicity values (Kalirai et al. 2008). 

To estimate the age of a given stellar population, the 
metallicity of the environment has to be taken into account. 
As it is well known, for a hxed initial mass, the lifetime of 
a star decreases with the metallicity, since less metallic - 
lower opacity - envelopes leads to higher luminosities. On 
the other hand, an increase on metallicity leads to an en¬ 
hanced mass loss during the giant phases of the evolution, 
resulting in lower hnal masses. In addition, other parame¬ 
ter characterizing white dwarf stars are expected to depend 
on metallicity. For instance, the amount of hydrogen left on 
a white dwarf star decreases for higher metallicity values 
(see e.g. Renedo et al. 2010). Also, a reduction of the initial 
metallicity tends to reduce the degree of chemical stratifica¬ 
tion in the central regions, which translates in a lower oxygen 
abundance at the center (Dominguez et al. 1999). This will 
affect the cooling rates, and in particular, the amount of en¬ 
ergy released by the crystallization process at Ten < 12 000 
K. 

In this work we focus on the study of the effects of dif¬ 


ferent metallicity values on the IFMR and on the lifetimes 
during the cooling sequences and its previous evolutionary 
stages. Specifically we center our analysis on the low and 
intermediate mass end of the IFMR, by considering masses 
ranging from ~ O. 8 M 0 to 2>Mq at the main sequence. To 
this end we compute full evolutionary sequences from the 
zero age main sequence (ZAMS), through the central stable 
hydrogen and helium central burning, the thermally puls¬ 
ing and mass-loss phases on the asymptotic giant branch 
(AGB), and Hnally to the white dwarf cooling phase. The 
computations were performed with the LPGODE evolution¬ 
ary code (Althaus et al. 2005a). We consider 7 metallicity 
values ranging from Z = 0.0001 to Z = 0.04, covering the 
values for halo and disk populations. 

This paper is organized as follows. Section [2] is devoted 
to present the input physics considered in the evolutionary 
computations. In section [3] we present our model grid and 
discuss the influence of metallicity on the pre-white dwarf 
evolution. We explore the influence of metallicity on the min¬ 
imum mass to reach helium burning at the beginning of the 
horizontal branch phase and the core mass growth during 
the TP-AGB phase. In section |4] we present and discuss our 
results for the theoretical IFMR for the different metallic¬ 
ity values and compare them with the observational data 
from open and globular clusters. Also, we compare our re¬ 
sults with other theoretical and semi-empirical determina¬ 
tions of the IFMR. In section [5] we analyze the influence 
of the metallicity on structural parameters of the resulting 
white dwarf models, such as the amount of hydrogen and 
the chemical composition of the core. Section | 6 ] is devoted 
to asses the effect of metallicity on the age determination 
of stellar populations. In particular we use our theoretical 
evolutionary sequences to study three well known globular 
cluster, NGG 6397, M4 and 47 Tuc. For this systems we es¬ 
timate the mass range at the cooling sequence. We present 
our concluding remarks in section 0 


2 NUMERICAL TOOLS 

The evolutionary computations presented in this work were 
calculated with an updated version of the LPGODE evolu¬ 
tionary code. The details on the code can be found in Al¬ 
thaus et al. (2005a, 2010), Renedo et al. (2010) and Romero 
et al. (2013). The LPGODE evolutionary code computes 
the complete evolution from the ZAMS, through the hy¬ 
drogen and helium burning stages, the thermally pulsating 
and mass loss stages on the asymptotic giant branch, to the 
white dwarf cooling evolution. This code has been used to 
produce very accurate white dwarf models (see Garcia-Berro 
et al. 2010; Althaus et al. 2010; Renedo et al. 2010 and ref¬ 
erence therein). The code has also been used to study the 
formation of extreme horizontal branch stars (Miller Berto- 
lami et al. 2008), the evolution of extremely low mass white 
dwarf stars (Althaus et al. 2013), and also the role of ther¬ 
mohaline mixing for the surface composition of low-mass red 
giants (Wachlin et al. 2011). 
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2.1 Pre-white dwarf evolution 

The LPCODE considers simultaneous treatment of non- 
instantaneous mixing and burning of elements (Althaus et 
al. 2003). The nuclear network accounts for 16 elements and 
34 reaction rates, that includes pp—chain, CNO—cycle, he¬ 
lium burning and carbon ignition (Althaus et al. 2005b). 
Reaction rates are taken from Caughlan & Fowler (1988) 
and Angulo et al. (1999). In particular, the ^^C(q, 7)^®0 
reaction rate, of special relevance for the carbon—oxygen 
stratification of the resulting white dwarf, was taken from 
Angulo et al. (1999). The radiative opacities are those of the 
OPAL project (Iglesias & Rogers 1996), including carbon— 
and oxygen—rich compositions, supplemented at low tem¬ 
peratures with the molecular opacities of Ferguson et al. 
(2005) and Weiss & Ferguson (2009). In LPCODF molecu¬ 
lar opacities are computed adopting the opacity tables with 
the correct abundances of the unenhanced metals and the 
appropriate carbon-oxygen ratio. For the present compu¬ 
tations, we have not considered carbon-enriched molecular 
opacities (Marigo 2002), which are expected to reduce ef¬ 
fective temperatures at the AGB (Wiess & Ferguson 2009). 
Conductive opacities are those from Cassisi et al. (2007). 
Neutrino emission rates for pair, photo and bremsstrahlung 
processes are taken from Itoh et al. (1996), and for plasma 
processes the treatment of Haft et al. (1994) is included. We 
adopted the standard mixing length theory with the free 
parameter a = 1.6. With this value, the present luminosity 
and effective temperature of the Sun, log Teg = 3.7641 and 
Lq = 3.842 X lO^^ergs”^ at an age of 4570 Myrs, are repro¬ 
duced by LPCODF when Z = 0.0164 and X — 0.714 are 
adopted, in agreement with the ZjX value of Grevesse & 
Sauval (1998), which are better reproduced by helioseismol¬ 
ogy than more recent determinations (Basu & Antia 2008). 

During the evolutionary stages prior to the ther¬ 
mally pulsing AGB phase, we allow the occurrence of 
extra—mixing episodes beyond each convective boundary 
following the prescription of Herwig et al. (1997). Note that 
the occurrence of extra—mixing episodes during the core he¬ 
lium burning phase largely determines the final chemical 
composition of the core of the resulting white dwarf star 
(Prada Moroni & Straniero 2002, Straniero et al. 2003). We 
treated the extra—mixing as a time-dependent diffusion pro¬ 
cess, assuming that the mixing velocities decay exponen¬ 
tially beyond each convective boundary. The diffusion co- 
efhcient is given by Dem = Do exp{—2z/fHp), where Hp 
is the pressure scale height at the convective boundary. Do 
is the diffusion coefhcient of unstable regions close to the 
convective boundary, and 2 is the geometric distance from 
the edge of the convective boundary (Herwig et al. 1997, 
2000). The free parameter / describes the efhciency of the 
extra—mixing process. It can take values as high as / 1.0, 

for over-adiabatic convective envelopes of DA white dwarfs 
(Freytag et. al 1996). However, for deep envelope and core 
convection / is expected to be considerably smaller because 
the ratio of the Brunt—Vaisala timescales of the stable to 
unstable layers decreases with depth. We assume / = 0.016, 
which accounts for the location of the upper envelope on 
the main sequence for a large sample of clusters and as¬ 
sociations (Schaller et al. 1992, Herwig et al. 1997, 2000). 
The same value accounts for the intershell abundances of 


hydrogen—deficient post—AGB remnants (see Herwig et al. 
1997; Herwig 2000; Mazzitelli et al. 1999). Finally, for the 
mass range considered in this work, the mass of the outer 
convection zone on the tip of the RGB only increases by 
1.5% if we change the parameter / by a factor of two. The 
breathing pulse instability occurring towards the end of the 
core helium burning is attributed to the adopted algorithm 
rather than to the physics of convection and therefore were 
suppressed in our computations (see Straniero et al. 2003 
for a detailed discussion). 

During the thermally pulsing AGB phase, 
extra—mixing episodes were disregarded. This may 
prevent the third dredge—up to occur in low mass stars. In 
particular, a strong reduction of extra—mixing episodes at 
the base of the pulse-driven convection zone seems to be 
supported by simulations of the s—process abundance pat¬ 
terns (Lugaro et al. 2003) and by observational inferences 
of the IFMR (Safaris et al. 2009). As a consequence, we 
expect that the mass of the hydrogen free core to gradually 
increase as evolution proceeds during this stage. Kalirai 
et al. (2014) estimated a growth of the core mass on the 
TP-AGB to be between 10% and 30% for initial masses be¬ 
tween I. 6 M 0 and 2 .OM 0 , that decreases steadily to ~ 10% 
at Mzams = 3.4Mq. However, third-dredge-up episodes 
are not completely suppressed in our models. In particular, 
the two more massive sequences for Z — 0.0001, those with 
initial masses 2.0 and 2.5 Mq do experience dredge-up 
episodes during the TP-AGB. We will comment on the 
effect of third dredge-up episodes on our determination of 
the final masses in section [331 In addition, the suppression 
of dredge-up episodes in our models prevents the formation 
of carbon rich stars, i.e, the G/O ratio remains below unity 
during all the TP-AGB stage. 

Mass loss was considered during red giant branch and 
core helium burning phases following Schroeder & Cuntz 
(2005). This mass—loss rate is based on the Reimers for¬ 
mula M = riLfR:t/Mt (Reimers 1975), but it contains 
two additional factors, one depending on the surface tem¬ 
perature and the other on the surface gravity of the star. 
This improves the agreement with the observed mass—loss 
rates for different types of stars, without the need to adjust 
the htting parameter rj (Shroeder & Cuntz 2005). We set 
rj = 8 X lO~'^'^M 0 yr~'^, that satisfies the well constrained 
RGB mass loss of globular cluster stars. 

During the AGB and TP-AGB phases, we employ the 
prescriptions of Groenewegen et al. (2009) and Vassiliadis & 
Wood (1993) for low and high metallicity sequences, respec¬ 
tively. The mass loss prescription presented by Vassiliadis 
& Wood (1993) for the TP-AGB phase is pulsation depen¬ 
dent, with an exponentially increasing rate for periods less 
than ~ 500 days and a constant superwind phase for periods 
beyond ~ 500 days, corresponding to a radiation—pressure 
driven wind. This prescription, or an updated version, is 
usually adopted to model the mass loss during the TP-AGB 
stage in evolutionary computations (e.g. Renedo et al. 2010, 
Karakas & Lattanzio 2007, Karakas 2010, Marigo & Girardi 
2007). In particular, Doherty et al. (2014) employed the 
mass—loss rate prescription of Vassiliadis & Wood (1993) 
to determine the IFMR of intermediate mass stars with ini¬ 
tial masses between 5.0 and 10.0 Mq for metallicities span- 


© 0000 RAS, MNRAS 000, 000-000 


4 A. D. Romero, F. Campos and S. O. Kepler 


ning the range of Z = 0.02 — 0.0001. Groenewegen et al. 
(2009) presented a mass loss prescription based on data of 
resolved stars on the Small and Large Magellanic Clouds 
obtained with Spitzer IRS that represents low metallicity 
populations. They found that, even when the Vassiliadis 
& Wood (1993) models described the data better, there 
are some deficiencies, in particular the maximum adopted 
mass—loss rate. The assumption of typical velocities and 
dust-to-gas ratios of galactic stars may lead to an overesti¬ 
mate of the mass—loss rate. Therefore we employ the mass 
loss prescription of Groenewegen et al. (2009) for low metal¬ 
licity sequences as it is better calibrated in this metallicity 
range. For metallicities larger than Z = 0.01, mass—loss 
rate from Groenewegen et al. (2009) are systematically lower 
than those from Vassiliadis & Wood (1993) and the well 
known superwind phase at the end of the thermally puls¬ 
ing stage is rarely reached even after the star moves sig¬ 
nificantly to redder effective temperatures. For this reason 
we employ the prescription of Groenewegen et al. (2009) for 
evolutionary sequences characterized by initial metallicities 
Z ^ 0.008; and the Vassiliadis & Wood (1993) mass—loss 
rates for sequences with initial metallicities Z ^ 0.01, that 
better reproduce mass loss rates for high Z. 


2.2 White dwarf evolution 

During the white dwarf evolution, we consider the distinct 
physical process that modify the chemical abundances dis¬ 
tribution. For instance, element diffusion strongly modifies 
the chemical composition profile throughout the outer lay¬ 
ers. Indeed, our sequences developed a pure hydrogen en¬ 
velope with increasing thickness as evolution proceeds. Our 
treatment of time—dependent diffusion is based on the mul¬ 
ticomponent gas treatment presented in Burgers (1969). We 
consider gravitational settling and thermal and chemical dif¬ 
fusion of H, ®He, ^He, ^^G, ^®G, and (Althaus et al. 
2003). To account for convection in our models we adopted 
the mixing length theory, in its ML2 flavor, with the free pa¬ 
rameter a = 1 (Tassoul et al. 1990). The metal mass fraction 
in the envelope is not fixed, but it is computed consistently 
according to the predictions of element diffusion. To account 
for the different envelope compositions, we considered radia¬ 
tive opacity tables from OPAL for arbitrary metallicities. 
For effective temperatures less than 10 000 K we include the 
effects of molecular opacity from the computations of Marigo 
& Aringer (2009). In addition, we considered the chemi¬ 
cal rehomogenization of the inner carbon—oxygen profile in¬ 
duced by Rayleigh—Taylor instabilities following Salaris et 
al. (1997). For effective temperatures below 10 000 K, outer 
boundary conditions are derived from non—gray model at¬ 
mospheres (Rohrmann et al. 2012). The equation of state is 
that of Segretain et al. (1994) for the high density regime, 
which accounts for all the important contributions for both 
the liquid and solid phases (Althaus et al. 2007), comple¬ 
mented at the low-density regime with an updated version 
of the equation of state of Magni & Mazitelli (1979). 

Finally, cool white dwarf stars are expected to undergo 
crystallization as a result of strong Goulomb interactions in 
their dense interiors (van Horn 1968). Grystallization occurs 
when the energy of the Goulomb interaction between neigh¬ 


boring ions is much larger than their thermal energy. This 
process leads to the release of latent heat and of gravitational 
energy associated with changes in the chemical composi¬ 
tion of carbon—oxygen profile (Garcfa—Berro et al. 1988ab). 
These additional energy sources associated with crystalliza¬ 
tion were included self-consistently, modifying the luminos¬ 
ity equation to account for both the local contribution of 
energy released from the core chemical redistribution and 
the latent heat, the later being of the order of 0.77fcsT per 
ion. To asses the enhancement of oxygen in the crystallized 
core we used the phase diagram presented in Horowitz et al. 
(2010) - see also Hugoto et al. (2012) and Schneider et al. 
(2012) -, consistent with the observations of Winget et al. 
(2009; 2010) of the white dwarf luminosity function of the 
globular cluster NGG 6397. Horowitz et al. (2010) employed 
molecular dynamic simulations considering liquid and solid 
phases in mixed composition, allowing a better determina¬ 
tion of the melting temperature and the composition of the 
liquid and solid phases. They predicted an azeotopic type 
phase diagram, and a melting temperature lower than that 
predicted by Segretain & Ghabrier (1993). Note that the re¬ 
sults obtained by Horowitz et al. (2010) are in good agree¬ 
ment with the results of Medin & Gumming (2010), based 
on a semi-analytic method. 

Althaus et al (2012) presented an study of the influence 
of the Horowitz et al. (2010) phase diagram on the evolution¬ 
ary properties of white dwarf stars, specifically on the cool¬ 
ing times. They found that the amount of mass redistributed 
in the phase separation is smaller for Horowitz et al. (2010) 
phase diagram than for the Segretain & Ghabrier (1993) 
one, leading to a smaller energy release from carbon—oxygen 
differentiation. Also, the composition changes are less sen¬ 
sitive to the initial chemical profile, so the cooling delay 
is less affected by the uncertainties in the carbon—oxygen 
initial abundances. Romero et al. (2013) study the effects 
of the two crystallization phase diagrams on the pulsation 
spectrum of variable DA white dwarf stars with stellar mass 
O.72M0. They showed that, for a given stellar mass, com¬ 
putations employing the Segretain & Ghabrier (1993) phase 
diagram was ~ 1000 K higher than that predicted by the 
Horowitz et al. (2010) formulation. From their asteroseis- 
mological fits they found that the phase diagram presented 
in Horowitz et al. (2010) is the one that better represents the 
crystallization process in the dense interiors of white dwarf 
stars. 

The LPGODE has been tested against other white 
dwarf evolution code, and the uncertainties in the white 
dwarf cooling ages arising from different numerical imple¬ 
mentations of stellar evolution equations were found to be 
below 2% (Salaris et al. 2013). 

3 MODEL GRID AND THE PRE-WHITE 

DWARF EVOLUTION 

3.1 Model Grid 

The current grid of evolutionary models is composed by a 
set of sequences with initial mass between O. 8 M 0 and 3 M 0 
at the ZAMS, considering also different metallicity values 
for the progenitor star between Z = 0.0001 to 0.04. Note 
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Figure 1. Hertzprung—Russell diagram of our evolutionary se¬ 
quences for Z = 0.004. From bottom to top: evolution of 
the O.92M0, 0.95Mq, 0.98Mq, I.OMq, 1.2M0, 1.5M0, 2.0Mq, 
2.25Mq, 2.5Mq and 3 .OM 0 model stars. The dashed line repre¬ 
sents the location of the Zero Age Main Sequence. Note that the 
sequences with an initial mass of O.92M0 and O.95M0 experi¬ 
ence an hydrogen sub-flash at high effective temperatures, before 
entering the cooling sequence. 


that in this mass range, all sequences end their lives as 
carbon—oxygen core white dwarf stars. The initial helium 
content for the starting model at the ZAMS was provided 
by the relation Y = 0.245 -I- 2Z, were Z represents the initial 
metallicity. Our results are summarized in Table [T] where 
we list the initial and final masses for all metallicity values 
considered in this work. For each sequence we also list the 
progenitor age, i.e., the age from the ZAMS to the point of 
maximum effective temperature before entering the white 
dwarf cooling curve (tprog). The sequences with Z = 0.01 
and Z = 0.001, except for the sequences with initial mass 
of 0.92 and I.IM 0 for the latter Z, are based on those from 
Renedo et al. (2010). Since these authors employ the Segre- 
tain et al. (1993) formulation of the phase diagram for crys¬ 
tallization, we recalculated the final evolution of the white 
dwarf cooling considering the Horowitz et al. (2010) pre¬ 
scription. The remaining sequences were computed specif¬ 
ically for this work. Additional sequences characterized by 
high and very low metallicity values, those with Z — 0.05 
and Z — 0.00001, were computed to asses the dependence of 
the final mass with metallicity, but also of other parameters 
as the amount of mass lost and the size of the helium core 
at the beginning of the central helium burning stage (see 
section 

The evolution in the Hertzprung—Russell diagram frm 
the ZAMS to advanced stages of white dwarf evolution for 
Z = 0.004 is shown in Figure [T] for sequences with initial 
masses larger than O.9OM0. Typically, sequence with initial 
mass lower than 0.95Mq experience an hydrogen sub-flash 
before entering its final cooling track. Hydrogen sub-flashes 



Figure 2. Hertzprung—Russell diagram of our evolutionary se¬ 
quences with initial stellar mass of IM 0 with metallicities of 
Z = 0.00001,0.0001,0.004,0.01,0.02 and 0.04 from top to bot¬ 
tom. 

are a consequence of an increase in the luminosity due to 
shell hydrogen burning at high effective temperatures dur¬ 
ing the evolution towards the white dwarf cooling curve. For 
Z — 0.004 the hydrogen sub-flashes are present for sequences 
with 0.9 and O.92M0. In addition, sequences characterized 
by the lowest initial stellar mass on the ZAMS usually do 
not go through the AGB and thermally pulsing AGB stages 
(O. 9 OM 0 for Z = 0.004). These sequences spend the core he¬ 
lium burning stage at the blue edge of the Horizontal Branch 
phase and move directly to the white dwarf cooling curve 
without experiencing the high luminosity phase. 

In Fig. [ 2 ] we show the evolutionary sequences charac¬ 
terized by a IM 0 at the ZAMS and different initial metal¬ 
licity values. When metallicity increases, the position of the 
main sequences moves towards lower luminosities and effec¬ 
tive temperatures, since the atmospheres of more metallic 
stars are more opaque, and then the larger opacities make 
the models move towards redder temperatures. On the other 
hand less metallic models are more luminous because the at¬ 
mospheres are more transparent. This is also true when we 
inspect the AGB and TP-AGB stages; less metallic models 
are dislocated to the blue because of the lower opacities. 
As we mention in section \m the mass-loss rate is metal¬ 
licity dependent (Vasiliadis & Wood 1993; Groenewegen et 
al. 2009), then the less metallic sequences have less efficient 
mass loss episodes, leading to a larger number of thermal 
pulses. 

In the present work, the hydrogen content of the se¬ 
quences correspond to the value predicted by standard stel¬ 
lar evolution and should be considered as upper limits for 
the maximum hydrogen content left in the white dwarf re¬ 
sulting for the evolution of single star progenitors. As we will 
show in section [5T] the amount of hydrogen left on a white 
dwarf envelope depends on stellar mass and metallicity (see 
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Table 1. Initial and final stellar masses (in solar units) and progenitor ages (in Gyr) for the metallicity values considered in Figure |4] 
Values corresponding to Z=0.01 and most to Z=0.001 are extracted from Renedo et al. (2010) (see text for details). 


Mzams 

Z = 

Mwd 

0.0001 

tprog 

Z = 

Mwd 

0.001 

tprog 

Z = 

Mwd 

0.004 

tprog 

Z = 

Mwd 

0.008 

tprog 

Z = 

Mwd 

0.01 

tprog 

Z = 

Mwd 

0.02 

tprog 

Z = 

Mwd 

0.04 

tprog 

0.80 

0.519 

12.742 













0.85 

0.534 

10.327 

0.505 

13.322 











0.90 

0.550 

8.421 



0.503 

11.898 

0.489 

15.767 







0.92 



0.536 

9.789 

0.515 

10.978 









0.94 







0.508 

13.037 







0.95 

0.561 

6.998 



0.524 

9.772 



0.493 

13.350 

0.488 

15.160 



0.98 





0.533 

8.736 

0.516 

11.679 







1.00 

0.569 

5.871 

0.553 

7.442 

0.537 

8.076 

0.524 

10.788 

0.525 

11.117 

0.511 

12.512 

0.511 

13.327 

1.05 







0.531 

8.053 







1.10 



0.566 

5.274 

0.547 

5.800 

0.538 

7.339 







1.20 





0.553 

4.310 

0.552 

5.047 







1.25 

0.621 

2.826 

0.593 

3.449 











1.50 

0.669 

1.593 

0.627 

1.897 

0.599 

2.080 

0.588 

2.382 

0.570 

2.700 





1.75 

0.708 

1.038 

0.660 

1.324 





0.593 

1.699 





2.00 

0.737 

0.752 

0.693 

0.947 

0.665 

1.074 

0.638 

1.074 

0.609 

1.211 

0.591 

1.286 

0.566 

1.368 

2.25 





0.691 

0.788 



0.632 

0.989 





2.50 

0.826 

0.423 



0.730 

0.579 



0.660 

0.742 





3.00 

0.875 

0.274 

0.864 

0.326 

0.817 

0.352 



0.705 

0.443 



0.656* 

0.461 


*Mass of the degenerate core at the last thermal pulse. 


e.g. Renedo et al. 2010, Romero et al. 2012). The amount of 
hydrogen left in a white dwarf star is an important factor 
to determine the cooling ages, since hydrogen has a larger 
opacity than helium. Cooling ages can be reduced by 10% for 
thin hydrogen envelopes (10“^°M*) at log(L/L 0 ) = —5.5 
(Catalan et al. 2008a) in comparison with thick ones. In 
addition, a thinner hydrogen envelope can be mixed with 
the underlying helium layer when the outer convective zone 
at the white dwarf atmosphere reaches the hydrogen-helium 
interface at low effective temperatures, turning the hydro¬ 
gen atmosphere star into a helium atmosphere white dwarf 
(Tremblay & Bergeron 2008; Romero et al. 2013; Kurtz et 
al. 2013). Evidence from asteroseismology (Castanheira & 
Kepler 2009; Romero et al. 2012, 2013) showed that the 
mass of the hydrogen layer spreads over 10 “^ — 10 “^°M« 
with a mean value around 10 “®M«, two orders of magni¬ 
tude smaller than the white dwarf models in the literature 
employed in evolutionary studies of this type of star. 


3.2 The core mass at helium ignition 

In this section we analyze the impact of metallicity on the 
mass of the helium core at the beginning of the central he¬ 
lium burning stage, and also on the total mass lost in the 
different giant phases. In Table [2] we list the final mass val¬ 
ues along with some other relevant parameters for all se¬ 
quences with initial mass of IMq and 2M©. We computed 
additional sequences, apart from those listed in Table[T]with 
Z = 0.00001, ^ = 0.03 and Z = 0.05 for Mzams = IMq, 
and Z = 0.03 for Mzams = 2Mq. The second column of 
Table m lists the mass of the helium core before the onset 
of helium burning, i.e., the hydrogen free core mass at the 
end of the RGB phase. For all sequences with a IM© pro¬ 
genitor and the sequences with a progenitor star with 2M© 


and metallicity Z ^ 0.01, this mass is the mass of the he¬ 
lium core at the maximum of the helium luminosity during 
the helium flash. Note that for sequences with initial mass 
IM©, the mass of the helium core increases with decreasing 
metallicity approximately 0.052M© in the range considered 
in this work. This can be seen also from Fig. [3] were we 
plot the mass of the helium core at the beginning of the 
central helium burning stage as a function of metallicity for 
sequences with initial stellar mass of IM© (squares) and 
2M© (circles). 

For sequences with progenitor mass of 2M© we have 
two groups, those sequences with initial metallicity below 
Z = 0.01 do not experience an helium ignition on a degen¬ 
erate core. As a consequence the mass of the helium core 
at the beginning of the central helium burning is signifi¬ 
cantly lower, around 0.32 — 0.33M©. On the other hand, 
sequences with initial metallicities ranging from 0.01 to 0.04 
ignite the central helium in a degenerate state. Although 
the mass of the helium core at this point is lower than that 
for the IM© sequences, it is higher compared to the metal 
poor cases. This change in the central helium core mass is 
clear in Fig. |3] In particular, we note that the helium core 
mass is maximum for Z = 0.02 and then decreases for lower 
and higher metallicities. Also, the maximum helium lumi¬ 
nosity at the moment of the helium flash is 3-4 orders of 
magnitude larger for the sequences with Z = 0.02. From 
stellar evolutionary theory we know that the size of the fi¬ 
nal hydrogen free core is mainly controlled by two process: 
mass-loss and overshooting. Mass-loss rates decreases with 
decreasing metallicity, allowing the core to increase its size 
given the larger amount of burning material. On the other 
hand, even though mass-loss rates increases with metallic¬ 
ity, the influence of overshooting is more important for high 
Z, also allowing the core to increase its size. For Z as high 
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Table 2. Different parameters as a function of the metallicity of the progenitor star with initial mass IMq (squares) and 2Mq (circles) 
at the ZAMS. The columns list the metallicity value, the mass of the hydrogen free core at the helium flash or at the beginning of the 
helium central burning stage, the mass lost during the Red Giant Branch and Asymptotic Giant Branch phases and the final mass at 
the white dwarf cooling sequence. All masses are in solar units. 


z 

Mhfc 

Mzams 

AM(RGB) 

= IMq 
AM(AGB) 

Mwd 

Mhfc 

Mzams 

AM(RGB) 

— 2iW0 
AM(AGB) 

Mwd 

0.00001 

0.5066 

0.1068 

0.3280 

0.564 





0.0001 

0.4986 

0.1374 

0.2960 

0.562 

0.3263 

0.0001 

1.2006 

0.738 

0.001 

0.4855 

0.1782 

0.2535 

0.553 

0.3227 

0.0021 

1.2725 

0.693 

0.004 

0.4807 

0.2342 

0.2168 

0.537 

0.3185 

0.0023 

1.3195 

0.665 

0.008 

0.4789 

0.2834 

0.1818 

0.524 

0.3225 

0.0028 

1.3492 

0.638 

0.01 

0.4789 

0.3218 

0.1410 

0.525 

0.4047 

0.0172 

1.3565 

0.609 

0.02 

0.4705 

0.3628 

0.1169 

0.511 

0.4172 

0.0243 

1.3777 

0.591 

0.03 

0.4661 

0.3870 

0.0883 

0.511 





0.04 

0.4596 

0.3912 

0.0887 

0.511 

0.3910 

0.0149 

1.4126 

0.566 

0.05 

0.4547 

0.3747 

0.1004 

0.515 







Figure 3. Mass of the hydrogen free core at the beginning of 
the central helium burning stage as a function of metallicity, for 
sequences with IMq (squares) and 2Mq (circles) at the ZAMS. 
Note that the larger the mass at the ZAMS, the smaller the mass 
of the hydrogen free core. 

as 0.05 the overshooting influence could be overcoming the 
influence of mass loss. 

Also listed in Table [2] are the total mass lost during the 
RGB and AGB phases. For the sequences with a progenitor 
star of IMq we see that the mass lost during the RGB in¬ 
creases with increasing metallicity, leaving a smaller amount 
of mass to be lost during the AGB phase. In fact, sequences 
with the highest values of Z do not experience any thermal 
pulses and the growth in the carbon-oxygen core is only 
due to the evolution during the AGB. On the other hand, 
the low metallicity sequences experience up to 4-5 thermal 
pulses, extending their stay at the high luminosity stage. For 
the sequences with a progenitor star of 2Mq the mass lost 
during the RGB is negligible compared with the total mass 
lost during the AGB and thermal pulsing AGB. Practically 


all mass is lost during the late stages of the evolution. Be¬ 
cause of this, the modeling of the mass-loss rates becomes 
so important during the AGB and TP-AGB. 

Finally, the dependence of the final mass with the pro¬ 
genitor metallicity also has an impact on the inner com¬ 
position of a given white dwarf star. Since the final mass is 
higher for lower metallicities, low mass carbon-oxygen white 
dwarf can have lower mass progenitors at lower metallicity. 
On the other hand, for lower main sequence stars, the mass 
of the helium core at the helium flash increases with de¬ 
creasing metallicity. For instance, a white dwarf with a stel¬ 
lar mass of ~ 0.5Mq resulting from a IMq progenitor star 
could have an helium core if the initial metallicity is lower 
than Z — 0.001, but it will be a carbon-oxygen core white 
dwarfs for higher metallicity values, according to our single 
star evolution computations. 

3.3 Core mass growth during the thermal pulsing 
AGB 

After helium exhaustion in the central core, the star moves 
to the AGB. In the early stages, the helium burning shell 
increases the carbon-oxygen core mass as it moves towards 
the surface. In our computations, overshooting is taken into 
account during the early AGB, with an efficiency param¬ 
eter of / = 0.016, while it is disregarded during the fol¬ 
lowing thermal pulsing stage. A reduction of extra—mixing 
episodes at the base of the driving convection zone during 
the TP-AGB could prevent the third dredge-up to occur, 
specially in low mass stars. The occurrence of the third 
dredge-up cause an instantaneous reduction of the core 
mass, dragging processed material, specially carbon, to the 
surface (see e.g. Wiess & Ferguson 2009; Karakas et al. 2010; 
Marigo et al. 2013). The efficiency of the third dredge-up is 
usually parametrized by the parameter iQ, which increases 
with each thermal pulse until a maximum value Amax (see 

^ The efficiency parameter A is defined as the fraction of the core¬ 
mass growth over the interpulse period that is dredge—up to the 
surface at the next thermal pulse 
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e.g. Marigo et al. 2013). The value of Amax typically in¬ 
creases with the stellar mass, while it decreases at larger 
metallicities. Then, dredge-up episodes are expected to be 
more efficient for high masses and low metallicities. Unfor¬ 
tunately, the efficiency A is one of the most uncertain pa¬ 
rameters of TP-AGB modeling, and it is very dependent on 
the adopted treatment of convection, mixing, and numerical 
codes (Marigo et al. 2012; Kalirai et al. 2014). 

Extra—mixing episodes during the TP-AGB also leads 
to the formation of carbon rich stars. Studies of AGB stel¬ 
lar populations in our Galaxy and the Magellanic Clouds 
shows that stars with initial mass < 2Mq and low metal- 
licity values do become carbon rich stars (Groenewegen et 
al. 1993; Stancliffe et al. 2005; Marigo & Girardi 2007). For 
instance, carbon rich stars must be formed for initial masses 
as low as IMq in the Large Magellanic Cloud in order to 
reproduce the carbon stars luminosity function (Sancliffe et 
al. 2005). Carbon enrichment of the atmosphere will cause 
an abrupt change in the molecular equilibrium leading to 
a rise of the atmospheric opacity. This causes an increase 
in the mass-loss rate, leading to an early end of the TP- 
AGB and thus a reduction of the final mass. Because in our 
computations we disregarded extra—mixing episodes, our se¬ 
quences in general do not experience dredge-up episodes 
during the TP-AGB evolution, except for the three more 
massive sequences, corresponding to the lowest metallicity 
of our grid. In our models, sequences with initial masses be¬ 
tween 2Mq and 3Mq and Z — 0.0001 do experience dredge- 
up episodes during the TP-AGB evolution, with values of 
Amax of ~ 0.22 — 0.36, increasing with stellar mass. The 
mass of the core for the onset of the third dredge-up also in¬ 
creases from 0.652M© to O.836M0 for the same mass range. 
This values are ~ 10 — 20% larger in comparison to the re¬ 
sults of Wiess & Ferguson (2009) for Z = 0.0005, adopting 
/ = 0.016. Finally, all our sequences that experience dredge- 
up episodes become carbon-rich stars. We found that it only 
takes one dredge-up episode to get a carbon-rich surface, 
leading to a C/0 ratio of ~ 8 — 10 at the end of the TP- 
AGB. This values of C/O are much higher than the C/O 
~ 1.28 — 2.0 found in Galactic carbon stars (Lambert et al. 
1986; Ohnaka et al. 2000). Note that the input physics in 
LPCODE is dedicated to model in detail the inner struc¬ 
ture of the star for all evolutionary stages, and specially the 
white dwarf stage. However, it does not compute a detailed 
atmosphere model, similar to most full evolutionary codes in 
the literature. It is no tailored to reproduce the atmosphere 
abundances found in AGB stars. 

To asses the effect of the extra—mixing episodes on our 
determinations of the Hnal mass, we calculated additional 
sequences allowing extra—mixing episodes to occur, by con¬ 
sidering different values of the overshooting parameter /. 
Wiess & Ferguson (2009) consider / = 0.016, and found 
that the hnal mass was very similar to that of the hrst 
thermal pulse, as a consequence of the enhanced dredge-up. 
On the other hand Lugaro et al. (2003) argued for a lower 
efficiency / = 0.008, to achieve a better agreement with 
detailed s—process observed abundance patterns. We com¬ 
puted sequences with initial mass of 2M© and metallicity 
Z = 0.01 with different values of / =0.016, 0.008 and 0.002. 
Note that the core mass growth during the thermal pulses is 


by far larger for sequences with initial masses ~ 2Mq , where 
the core mass at the Hrst thermal pulse shows a strong min¬ 
imum (Renedo et al. 2010; Marigo et al. 2013; Kalirai et 
al. 2014). All three sequences computed with some amount 
of overshooting do become carbon stars after 10 — 11 ther¬ 
mal pulses and ~ 4 dredge-up episodes, with G/O ratios of 
~ 1.3 — 1.5 at the end of the TP-AGB evolution, and val¬ 
ues of Amax from 0.63 for / = 0.002 to 0.94 for / = 0.016. 
Regarding the core mass growth, we found that the Hnal 
mass of the hydrogen free core is reduced by 6.2%, 7.6% and 
9.4% when / takes values of 0.002, 0.008 and 0.016, respec¬ 
tively. As a result, the Hnal mass can be overestimated by 
less than 10% in our computations for this metallicity. For 
lower values of Z, like Z — 0.004 and the same initial mass, 
we found a larger reduction of the Hnal mass by ~ 15%, when 
/ = 0.008. This is in agreement with the fact that dredge- 
up episodes are more efficient for lower metallicity values. 
Karakas et al. (2002) showed that even for low metallicity 
values as % = 0.004, Amax is maximum for Minitiai ~ 2.5Mq 
and then slowly decreases for larger masses, leading to a 
less efficient reduction of the core mass due to dredge-up 
episodes. Then, in the worst case scenario, the uncertainties 
in Hnal masses from our treatment of overshooting episodes 
might be up to ~ 15% for lower metallicities. Note that, ad¬ 
ditional uncertainties on the Hnal core mass arise from the 
possible reduction of TP-AGB lifetimes when C-rich molec¬ 
ular opacities are considered, leading to an increase in the 
surface opacity and an enhanced mass-loss rate at this stage. 

Finally, Kalirai et al. (2014) performed an study of the 
core mass growth covering the intermediate-mass range of 
the IFMR. They use data from four open clusters with nearly 
solar metallicities: NGG 6819, NGG 7789, Hyades and Prae- 
sepe, the Hrst three included in our analysis of the initial- 
to-Hnal mass relation (see section|4ll. As a result, they found 
that, in order to Ht the data, the dredge-up efficiency Amax 
increases from zero to 0.5 for initial mass from ~ 2M© 
to ~ 3 Mq and then progressively decreases for larger ini¬ 
tial masses. From our theoretical results with Z = 0.002 
we found a core mass growth of 10% to ~ 28% for stars 
with initial mass 1.6Mq and 2.8Mq, decreasing steadily to 
~ 6% at Minitiai = 3.7M0, somewhat smaller than the re¬ 
sults of Kalirai et al. (2014), but still in excellent agreement. 
In addition, Kalirai et al. (2014) found that the dominant 
process governing the growth of the core is largely the stel¬ 
lar wind, while the third dredge-up, although not negligible, 
plays a secondary role. They also found that the mass loss 
prescription that better Hts the semi-empirical data is that 
of Vasiliadis & Wood (1993), the same mass loss prescription 
employed in our work for solar metallicity. 


4 INITIAL-TO-FINAL MASS RELATION 

In this section we show our results for the theoretical IFMR 
and its dependence with the diflerent metallicity. Given our 
model grid, we focus on the low and intermediate mass do¬ 
main of the IFMR, and consider sequences with an initial 
mass ^ 3Mq. This accounts for globular clusters and old 
open clusters, having main sequence stars within this mass 
range. For young open clusters with ages of the order of 
~ 100 Myr, low mass stars have no time to evolve out of 
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Figure 4. Theoretical IFMR for different metallicity values from Z = 0.0001 to Z = 0.04 and initial mass ^ 3Mq. Values for Z = 0.01 
and initial mass larger than 3Mq are taken from Renedo et al. (2010). All computations were made by using the LPCODE evolutionary 
code. We also include some observational data for globular and old open clusters — see text for details. (A color version of this figure is 
available in the on-line journal.) 


the main sequences yet. The initial mass (Mzams), the final 
mass (Mwd) and the progenitor lifetime for each evolution¬ 
ary sequences for all the metallicities are listed in Table [T] 

We consider the final mass as that of the white dwarf 
star in the cooling curve, where the mass-loss rate becomes 
negligible, not the mass of the hydrogen free core in the first 
thermal pulse. As show in section 13.31 and by Kalirai et al. 
(2014) the mass of the hydrogen free core increases during 
the TP-AGB as much as 30% for Minitiai ~ 2.0Mq. Then, 
the IFMR for a fixed metallicity by the end of the thermally 
pulsing AGB becomes markedly different from that deter¬ 
mined by the mass of the hydrogen free core at the first 
thermal pulse (Althaus et al. 2010). For instance, if we take 
a white dwarf with a stellar mass of O.63M0, the initial mass 
would be overestimated by ~ 0.75Mq if we were to consider 
an IFMR at the first thermal-pulse, instead of the real one, 
at advanced stages in the thermally pulsing AGB. For the 
carbon-oxygen composition expected in the core of a white 
dwarf, this is an important issue (Althaus et al. 2010). 

In Figure |4] we plot the initial-to-final mass relation. 
Each curve represents a set of theoretical computations with 
a given initial metallicity value. From this figure we can see 


that the IFMR shows a dependence with metallicity: for a 
given initial mass of the progenitor star, the mass of the re¬ 
sulting white dwarf remnant decreases when metallicity in¬ 
creases. This effect is noticeable for progenitor mass above 
~ IMq. For instance, the mass of the resulting white dwarf 
evolving from a 2Mq progenitor star is ~ 21% larger when 
the metallicity decreases from Z — 0.0001 to 0.01. We ex¬ 
pect the difference to be even larger for more massive stars. 
For low mass progenitor the differences are smaller, but still 
not-negligible. For a star with an initial mass at O.S5Mq 
there is an increase in the white dwarf mass of ~ O.O3M0 
when the metallicity decreases from 0.001 to 0.0001. Even 
if the differences in the final mass are not very large, the 
differences in the time spent during the stages previous to 
the white dwarf sequence will be, because the evolutionary 
time scales for low mass stars are very long. This age differ¬ 
ence is important when we use evolutionary computations 
to estimate the age of stellar populations. 

For comparison purposes, we include in Figure [T] ob¬ 
servational data corresponding to open and globular clus¬ 
ters with a Main Sequence Turn Of mass in the range 
~ 0.8 — 4.OM0. We consider the data from two globular 
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Table 3. Results from our fittings for the initial and final masses for the white dwarf stars corresponding to five different clusters. 
Effective temperature and superficial gravity values for NGC 6819, NGC 7789 and Hyades are taken from Kalirai et al. (2014), for NGG 
6791 are taken from Kalirai et al (2007), for NGC 6397 are taken from Moehler et al. (2004) and for M4 are from Kalirai et al. (2009). 


Cluster 

star 

Tea (K) 

logs 

Mwd/Mq 

teool{Myrs) 


NGG 6397 

WF4-358 

18800 ± 340 

7.72 

± 

0.06 

0.511 

± 

0.023 

75.75 ± 8 

0.829 

± 

0.030 

M4 

WD-00 

20600 ± 600 

7.75 

± 

0.09 

0.522 

± 

0.027 

52.89 ± 8 

0.885 

± 

0.011 

M4 

WD-04 

24300 ± 500 

7.68 

± 

0.08 

0.516 

± 

0.028 

27.48 ± 3 

0.885 

± 

0.011 

M4 

WD-06 

25600 ± 500 

7.87 

± 

0.08 

0.572 

± 

0.036 

21.08 ±5 

0.885 

± 

0.011 

M4 

WD-15 

24300 ± 600 

7.79 

± 

0.09 

0.534 

± 

0.039 

27.29 ±3 

0.885 

± 

0.011 

M4 

WD-20 

19700 ± 600 

7.73 

± 

0.10 

0.519 

± 

0.038 

58.83 ± 17 

0.885 

± 

0.011 

M4 

WD-24 

25700 ± 500 

7.70 

± 

0.07 

0.522 

± 

0.018 

22.32 ±3 

0.885 

± 

0.011 

NGG 7789 

NGC 7798-5 

31600 ± 200 

7.89 

± 

0.05 

0.636 

± 

0.024 

8.66 ±0.3 

1.908 

± 

0.051 

NGC 7789 

NGC 7798-8 

25000 ± 400 

8.06 

± 

0.07 

0.663 

± 

0.037 

26.15 ±9 

1.917 

± 

0.051 

NGC 6819 

NGC 6816-6 

21900 ± 300 

7.89 

± 

0.04 

0.532 

± 

0.017 

43.15 ±4 

1.561 

± 

0.050 

NGC 6819 

NGC 6816-7 

16600 ± 200 

7.97 

± 

0.04 

0.599 

± 

0.010 

141.63 ± 15 

1.585 

± 

0.054 

Hyades 

WD0325±096 

14670 ± 380 

8.30 

± 

0.05 

0.786 

± 

0.030 

357.59 ± 69 

3.655 

± 

0.692 

Hyades 

WD0406±169 

15810 ± 290 

8.38 

± 

0.05 

0.836 

± 

0.019 

337.26 ± 43 

3.682 

± 

0.536 

Hyades 

WD0421-I-162 

20010 ± 320 

8.13 

± 

0.05 

0.692 

± 

0.028 

97.75 ± 23 

2.906 

± 

0.109 

Hyades 

WD0425±168 

25130 ± 380 

8.12 

± 

0.05 

0.696 

± 

0.027 

30.48 ± 6 

2.804 

± 

0.109 

Hyades 

WD0431±126 

21890 ± 350 

8.11 

± 

0.05 

0.684 

± 

0.028 

61.76 ± 13 

2.851 

± 

0.096 

Hyades 

WD0437±138 

15120 ± 360 

8.25 

± 

0.09 

0.757 

± 

0.055 

305.00 ± 59 

3.498 

± 

0.377 

Hyades 

WD0438±108 

27540 ± 400 

8.15 

± 

0.05 

0.717 

± 

0.028 

17.65 ±5 

2.785 

± 

0.083 

Hyades 

WD0348±350 

14820 ± 350 

8.31 

± 

0.05 

0.792 

± 

0.029 

363.95 ± 73 

3.836 

± 

0.709 

Hyades 

HS0400±1451 

14620 ± 60 

8.25 

± 

0.01 

0.756 

± 

0.068 

325.17 ± 23 

3.612 

± 

0.403 

Hyades 

WD0625±415 

17610 ± 280 

8.07 

± 

0.05 

0.654 

± 

0.027 

137.70 ±21 

2.967 

± 

0.106 

Hyades 

WD0637±477 

14650 ± 590 

8.30 

± 

0.06 

0.786 

± 

0.036 

361.94 ± 111 

3.824 

± 

0.923 

NGC 6791 

WD7 

14800 ± 300 

7.91 

± 

0.06 

0.565 

± 

0.030 

196.66 ± 40 

1.166 

± 

0.028 


clusters, M4 and NGC 6397, and four open clusters, NGC 
6791, NGC 6819, NGC 7789 and Hyades. The cluster and 
star identifications and the spectroscopic parameters for the 
star sample are listed in the first four columns of TableO Ef¬ 
fective temperature and superficial gravity values were taken 
from the works of Kalirai et al. (2014) for NGC 6819, NGC 
7789 and Hyades, from Kalirai et al. (2007) for NGC 6791, 
from Moehler et al. (2004) for NGC 6397 and from Kalirai et 
al. (2009) for M4. These spectroscopic parameters were de¬ 
rived from the well know technique which involves fitting the 
Balmer lines shape of the spectrum to model atmospheres 
(e.g. Bergeron et al. 1992). With Teu and logy constrained, 
we estimated the spectroscopic mass values (column 5 of 
Table [3| by a linear interpolation of carbon-oxygen core 
cooling tracks in the log g — Tea, computed from full evolu¬ 
tionary models. The cooling ages (tcooi) for the white dwarf 
stars were extracted also from our evolutionary tracks (col¬ 
umn 6 from Table [S]). By subtracting tcooi from the total 
age of the cluster, we obtained the lifetime of the progenitor 
star (tprog) at the point of highest temperature in the post- 
AGB stage, before entering the white dwarf cooling curve. 
We assumed that the total ages of the clusters are 2.5 ± 0.2 
Gyr for NGC 6819 and 1.4 ±0.1 Gyr for NGC 7789 (Kalirai 
et al. 2001, 2008), 625 ± 50 Myrs for the Hyades (Perryman 
et al. 1998; Claver et al. 2001), 8.5 ± 0.5 Gyr for NGC 6791 
(Kalirai et al. 2007), 11.7 ± 0.3 Gyr for NGC 6397 (Hansen 
et al. 2013) and 11.6 ± 0.6 Gyr for M4 (Bedin et al. 2009). 
The initial mass values (last column of Table |31) have been 
computed directly from the pre-white dwarf lifetime tprog, 
using the mass-lifetime relations from pre-white dwarf the¬ 
oretical evolutionary models. For each cluster we take into 


account the metallicity value and consider a grid of progeni¬ 
tor sequences with the closest Z to compute the initial mass 
from the progenitor age. 

Note that there is a good agreement between theoret¬ 
ical IFMR and the semi empirical values. In general, all 
the points representing the star sample listed in Table [3] 
are accounted for by the range covered by our theoretical 
curves, within the metallicity range considered by our com¬ 
putations. Then, the spread observed in the data can be 
solved by considering a spread in the metallicity values. A 
more detailed analysis of Figure [T] shows that for stars be¬ 
longing to M4 and NGC 6397, at the low mass limit, and 
to the Hyades, for high mass values, the theoretical curves 
with the right metallicity fits the points within the uncer¬ 
tainties. For NGC 6819, the data is better represented with 
somewhat lower metallicities than the listed observational 
value Z = 0.017. However, asteroseismic inferences from 
red giants stars from this cluster suggest that metallicity 
might be significantly lower than solar, between Z = 0.0033 
and Z = 0.008 (Hekker et al. 2011). For the very metal 
rich open cluster NGC 6791, represented by a single ob¬ 
ject WD7, and the near solar metallicity open cluster NGC 
7789, the accordance between the semi empirical data and 
the theoretical IFMR is poor. For NGC 7789 the metallic¬ 
ity value Z = 0.014 is an average of literature values ob¬ 
tained by means of different photometric and spectroscopic 
techniques (see Wu et al. 2007), ranging from Z = 0.0085 
{[Fe/H] = —0.35) to solar C Z = 0.019). Our results would 
be in better agreement with a metallicity value closest to 
the lower limit ~ 0.008. In the case of NGC 6791 we ex¬ 
pect a higher mass-loss rate given its high metallicity. How- 
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Figure 5. Theoretical Initial-to-Final mass relation for sequences 
with metallicity values Z = 0.004 (black circles), Z = 0.008 (red 
squares), Z = 0.02 (green triangles up) and Z = 0.04 blue trian¬ 
gles down). Hollow symbols correspond to the results of Weiss & 
Ferguson (2009) while fill symbols correspond to results obtained 
in this work - see text for details. We consider the initial mass 
range where the two studies overlap. (A color version of this figure 
is available in the on-line journal.) 


ever the final mass obtained for WD7 is around ~ O.OSMq 
higher than the expected value, indicating that the theo¬ 
retical mass-loss rates at high metallicity values should be 
smaller than assumed, in agreement with the results ob¬ 
tained by Miglio et al. (2012) for this cluster. In summary, 
different metallicty values of the progenitor stars leads to 
different final masses on the white dwarf colling curve, caus¬ 
ing the observed spreading. 


4.1 Comparison with theoretical and 
semi-empirical IFMR 

In this section we compare our results with other theoretical 
determinations of the IFMR performed using independent 
evolutionary codes. In particular our definition of the final 
mass is the real mass of the white dwarf star at the cooling 
curve and not the mass of the hydrogen free core at the first 
thermal pulse. As it was shown by Kalirai et al. (2014) the 
mass of the core increases during the TP-AGB, resulting in 
a final mass larger than the one at the first thermal pulse. 
Therefore we select theoretical computations that consider 
the AGB evolution in the definition of the IFMR to compere 
with our results. 

Weiss & Ferguson (2009) computed sequences from the 
ZAMS to the end of the TP-AGB for 5 values of metallicity 
between Z = 0.0005 and Z — 0.04 and initial masses ranging 
from 1.0 Mq to G.OMq. We compare our results with those of 
Weiss & Ferguson (2009) in Fig. Only the metallicity val¬ 



Figure 6. Theoretical Initial-to-Final mass relation computed 
in this work (gray area) compare to the semi-empirical relations 
from Weidemann (2000, full line), Catalan et al. (2008, dashed 
line), Salaris et al. (2009, dashed-dot-dot line) and Kalirai et al. 
(2008, dashed-dot line). 


ues and the initial mass range where both studies overlap are 
plotted. At the low mass end we have good agreement, spe¬ 
cially for high metallicity progenitors. Differences are more 
important for high mass stars and low metallicity progeni¬ 
tors. In their calculation, envelope overshooting was allowed 
during the TP-AGB, with / = 0.016 for all sequences. This 
overshooting prescription leads to an enhanced reduction of 
the final core mass, leading to an IFMR very close to the 
relation between the initial mass and the core mass at the 
first thermal pulse. Note the minimum in the final mass for 
Weiss & Ferguson (2009) around ~ 2Mq corresponding to 
the limit between the sequences that experience an helium 
flash at the beginning of the central helium burning phase, 
and those that do not (see section lT^ . This minimum is also 
present in our computations when we consider the mass of 
the hydrogen free core at the first thermal pulse as the final 
mass (see Renedo et al. 2010). 

Other authors (e.g. Karakas 2010; Marigo & Gi- 
rardi 2007; Marigo et al. 2013) have computed detailed 
atmosphere models of TP-AGB evolution, considering 
parametrized overshooting treatments and also a detailed 
composition of the surface. They calibrate the dredge-up on¬ 
set and efficiency with observations of carbon surface abun¬ 
dances of AGB stars. Because in Karakas (2010) and Marigo 
& Girardi (2007) the authors consider dredge-up episodes 
during the TP-AGB we expect final mass for sequences with 
initial masses above ~ 2 M 0 to be higher in our computa¬ 
tions. We showed in section 13.31 that for Z ~ 0.02 our final 
masses can be overestimated by less than 10 %, while for 
lower metallicity this difference can reach ~ 20 %. 
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Karakas (2010) computed a large grid of TP-AGB mod¬ 
els covering the metallicity range between Z = 0.0001 and 
Z = 0.02, and initial masses between 1 and 6Mq. She con¬ 
sidered a detailed treatment of the third dredge-up in order 
to reproduce the atmospheric carbon abundances observed 
in held stars. As we expected, for initial mass larger than 
~ 1.8 — 2.0Mq we obtain hnal masses larger than those 
of Karakas (2010), up to a ~ 12% for Z = 0.0001 and 
Afinitiai = 3M©. For low initial masses our models show sig- 
nihcantly lower hnal masses than that of Karakas (2010), up 
to ~ 27% for IM© and Z = 0.0001. This large difference can 
be mainly due to the lack of mass loss during the pre-AGB 
evolution adopted for the initial models. As in Wiess & Fer¬ 
guson (2009), we consider mass loss during the RGB phase 
modeled by the Reimers formulae (Reimers 1975). There¬ 
fore our sequences begin the AGB evolution with a smaller 
mass than that of the ZAMS, limitng also the number of 
thermal pulses and thus, the growth of the hydrogen free 
core during the TP-AGB. Similar results to those of Karakas 
(2010) are found when we compare our IFMRs with those 
of Marigo & Girardi (2007). They also computed detailed 
models for TP-AGB evolution for Z — 0.0001 to Z = 0.03 
and Minitiai = 0.5 — 5.0M©, calibrating the onset and the 
efficiency of the third dredge-up with observations, such as 
carbon star luminosity functions in the Magellanic Clouds 
and TP-AGB lifetimes in Magellanic star clusters. The ini¬ 
tial models at the first thermal pulse are taken from Girardi 
et al. (2000), derived again from constant mass computa¬ 
tions. They still find a significant dependence with metal¬ 
licity, leading to an allowed final mass range of O.IM© wide 
for a given initial mass (see their figure 26). 

Finally, Kalirai et al. (2014) found a grid of best fit 
models from the model grid of Marigo et al. (2013), therefore 
without mass loss in the RGB, for Z = 0.02 that better 
represents the core mass growth during the TP-AGB for 
solar metallicity in the initial mass range of 1 — 4.4M©. The 
adopted dredge-up efficiencies varies with the initial mass. 
From their best-fitting TP-AGB models we can extract an 
IFMR (see their Table 3) and compare their predictions with 
ours. As a result, we found that the final masses from Kalirai 
et al. (2014) are somewhat larger than ours in the whole 
initial mass range. The differences are small, up to 4% 
compared to our Z = 0.01 grid, and 8% as compared to 
our Z = 0.02 models. The smaller differences are for initial 
mass around 2.5M©, where dredge-up efficiencies are larger 
for this metallicity. 

In addition to theoretical computations, several studies 
where dedicated to determine a semi-empirical IFMR from 
observations of star belonging to globular and open clusters 
and binary systems on the Galactic field. We selected the 
semi-empirical IFMR from Gatalan et al. (2008b), Salaris et 
al. (2009) and Kalirai et al. (2008) to compare with our the¬ 
oretical results. We also include the classical semi-empirical 
IFMR from Weidemann (2000). The results are shown in 
Figure[6l Most of the semi-empirical results fall between our 
theoretical IFMR characterized by Z = 0.004 and Z = 0.02. 
At the low mass end, below ~ 1.7M©, the IFMR of Catalan 
et al. (2008b) and Kalirai et al. (2009) fall between the 
Z = 0.004 and Z = 0.01 curves. Salaris et al. (2009) IFMR 
is better matched by our Z = 0.02 curve. At the high mass 


end, the semi-empirical IFMR are better matched by con¬ 
sidering metallicity values between Z = 0.008 and Z = 0.01. 
This is a consequence of the sample used to determine the 
IFMR, corresponding to younger, and more metallic open 
clusters, and common proper motion pairs from the solar 
vicinity. In particular, Kalirai et al. (2008) extended the low 
mass end towards lower initial masses by including the old 
metal-rich open cluster NGC 6791. In all cases, the semi- 
empirical IFMR demonstrate that lower mass main sequence 
stars produce lower mass white dwarf stars (Kalirai 2013). 


5 WHITE DWARF MODELS AND THE 

METALLICITY DEPENDENCE 

5.1 Hydrogen envelope thickness 

As shown by Iben & MacDonald (1986) -see also Renedo 
et al. (2010); Miller Bertolami et al. (2013)- the amount of 
hydrogen left in a white dwarf star is metallicity dependent, 
being thicker for models with less metallic progenitors. As a 
result of the larger hydrogen envelopes, residual H burning 
is expected to become more relevant in white dwarf stars 
with low metallicity progenitors. Results from Miller Berto¬ 
lami (2013) showed that for metal poor stars (Z 0.0001), 
the luminosity is completely dominated by nuclear burning, 
even at rather low luminosities. This leads to a significant 
delay in the cooling times, as compared to stars with higher 
metallicities, such as Z = 0.01, in which nuclear burning 
does not play a leading role, and most of the energy release 
come from the thermal energy stored in the interior. From 
the top panel of Figure[7l where we show the total hydrogen 
content in solar units in terms of the white dwarf mass at 
the point of maximum effective temperature, we see that, 
in general, a lower hydrogen content is found for a higher 
metallicity value for the progenitor star. Some departures of 
the expected trend can be noticed from this figure also. For 
instance, the change in the slope for high mass Z = 0.0001 
models is due to the occurrence of the third dredge-up for 
the two more massive models, as also shown in Figure 1 
of Miller Bertolami (2013). For Z = 0.001, the hydrogen 
content for the 0.627M© model is lower than that of the 
Z=0.004. The largest departures from the general trend are 
found for low mass white dwarfs, in particular those with 
stellar masses below ~ 0.53M©. As we mention in sec. 13.11 
sequences with the lowest masses of the grid experience hy¬ 
drogen sub-flashes before entering its final cooling track. In 
particular, the sequence with initial metallicity Z = 0.02 and 
white dwarf mass 0.488M© experience a late thermal pulse 
which significantly reduce the amount of hydrogen left in 
the star (Althaus et al. 2005b). 

As residual hydrogen burning is active, the mass of the 
hydrogen content decreases, until it reaches a certain value 
(-^ 8 X 10“®M© for solar metallicity stars, Renedo et al. 
2010 ) where the pressure at the bottom of the envelope is 
not large enough to support further nuclear reactions, and 
the hydrogen content reaches a stationary value. In addition, 
time-dependent chemical diffusion builds a pure hydrogen 
envelope when the star reaches lower effective temperatures 
(~ 30 000 K for solar metallicity). Then, at late stages of 
the cooling, a smoother distribution of the hydrogen mass 
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Figure 7. Hydrogen mass at the cooling sequence as a function 
of the white dwarf mass, at the point of maximum effective tem¬ 
perature at the beginning of the white dwarf cooling branch (top 
panel) and for models at late stages on the white dwarf cooling 
curve, with effective temperature ~ 12 000 K. Each curve corre¬ 
spond to a different initial metallicity (see references on the plot). 
A color version of this figure is available in the on-line journal. 

in terms of the stellar mass and metallicity can be expected. 
This can be seen for our model grid from the bottom panel 
of Figure [T] where we depict the total hydrogen mass as a 
function of the white dwarf mass at ~ 12 000 K. In particu¬ 
lar, at the low mass end, the hydrogen envelope mass clearly 
decreases as the initial metallicity increases. 

A comparison between the resulting hydrogen content 
for a fixed white dwarf mass shows a spread with metal¬ 
licity. For instance, for a white dwarf with a stellar mass 
of ~ O. 57 OM 0 we found that the reduction of the hydro¬ 
gen content due to an increase in the initial metallicity from 
Z = 0.001 to Z = 0.04 can be up to a factor of 2. This differ¬ 
ence is small and do not introduce a measurable difference 
on the cooling time scale at low luminosities. However we 
must take into account that the hydrogen content resulting 
from canonical evolutionary computations is a higher limit. 
Asteroseismological studies show that the mean hydrogen 
mass for pulsating ZZ Cetis stars is ~ 10“®M*, around two 
orders of magnitude lower than the canonical values for most 
of the sequences computed in this work (Castanheira & Ke¬ 
pler 2009; Romero et al. 2012, 2013). 

5.2 Central oxygen abundances 

Figure [8] displays the central oxygen abundance by mass 
left after core helium burning in term of the white dwarf 
mass for three metallicity values Z — 0.0001, Z — 0.004 and 
Z = 0.01. The numbers in the figure indicate the mass of the 
progenitor star for some points. The predictions of our calcu¬ 
lations for Z — 0.01 are in good qualitative agreement with 



Figure 8. Central oxygen abundance in terms of the white dwarf 
mass, for three metallicities, Z = 0.0001 (black circles), Z = 
0.004 (red squares) and Z = 0.01 (blue triangles). The values are 
taken after rehomogeneization but before crystallization sets in 
the central core. The numbers indicate the progenitor mass at 
the ZAMS for some points. 

the computations of Salaris et al. (1997) for solar metallicity, 
as shown in Althaus et al. (2010). We recall that the final 
carbon—oxygen stratification of the emerging white dwarf 
depends on both the efficiency of the ^^C(a, 7)'^®0 reaction 
rate and the occurrence of extra—mixing episodes toward 
the late stages of core helium burning. Note that for all sets, 
a maximum in the oxygen abundance is predicted, being 
the mass of the white dwarf model larger for lower metallic¬ 
ity values. However, when we consider the initial mass, we 
found that the maximum is located approximately at the 
same value, around 2.25Mq for Z = 0.004 and Z = 0.01 
and around 2Mq for Z = 0.0001. This maximum is also 
present in Salaris et al. (1997) models. After that, the oxygen 
abundance decreases and tends to a stable value oscillating 
between 0.61 and 0.63, as also shown in Romero et al. (2013) 
for white dwarf models with Z = 0.01. Finally, note that the 
changes in the oxygen and carbon central abundances due 
to metallicity can get up to ~ 10% in the metallicity range 
shown in Fig. [S] 


6 THE AGE OF STELLAR POPULATIONS 
AND THE METALLICITY PARAMETER 

As it is known, the age of the model before entering the final 
cooling sequence is metallicity dependent; for a given pro¬ 
genitor mass on the ZAMS, the age of the pre-WD star at 
the point of maximum effective temperature is significantly 
smaller for lower metallicities. For progenitor stars corre¬ 
sponding to initial masses above 1.5 Mq, the evolutionary 
scales are small, from 2 Gyr to hundreds of Myr. But for 
low mass stars, those that present a fraction of its core in 
a degenerate state after central hydrogen exhaustion, the 
evolutionary scales are much longer, and can reach 15 Gyr. 
In addition, for low mass stars the time spent on the RGB 
phase is not negligible and can reach up to ~ 2 — 3 Gyr 
for IMq. Thus the metallicity dependence of the pre-WD 
age will be more important for low mass sequences and older 
populations like globular clusters. We can see from Table [T] 
that the pre-WD age of a IMq progenitor star is reduce from 
12.5 Gyr to 5.9 Gyr when the metallicity goes from a solar 
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Figure 9. Lifetime previous to the white dwarf stage versus stel¬ 
lar mass on the cooling curve. Each curve corresponds to a single 
initial metallicity. (A color version of the figure is available on the 
on-line version). 


value {Z = 0.02) to very metal poor values characteristic of 
halo population {Z = 0.0001). 

Figure O shows the pre-WD age (tprog) in terms of the 
white dwarf mass for all metallicity values considered in this 
work. The general trend shows that the age, for a given stel¬ 
lar mass, decreases with metallicity (e.g. Isern et al. 2005). 
When the step on metallicity is sufficiently small, we can 
notice a departure from the general trend, caused by the 
difference in the progenitor mass. As it was shown in section 
[4] the mass of the resulting white dwarf star decreases when 
metallicity increases (see Table [T]). This is a consequence of 
an enhanced mass loss rate with increasing metallicity dur¬ 
ing the giant stages. Note that the metallicity dependence 
of the final mass will still be present if we were to consider 
the mass of the hydrogen free core at the first thermal pulse, 
since the mass of the carbon-oxygen core after helium burn¬ 
ing increases with decreasing metallicity, as shown in Figure 
[31 Two effects are present having opposite influences. For a 
fixed white dwarf mass, the lifetime of the progenitor star 
decreases with metallicity but, on the other, the mass of the 
progenitor is smaller and its lifetime increases (e.g. Isern et 
al. 2005). For a low mass white dwarf, the age difference can 
be as high as 8 Gyr. 

In Figure [in] we show the age of the progenitor star 
in terms of the initial metallicity. From this figure the dif¬ 
ferences on the white dwarf total age due to metallicity is 
clearer. We found that the age decrease for low Z but the 
dependence is not always monotonous, indicating the influ¬ 
ence of another parameter, i.e. the IFMR. As the stellar 
mass increases, the pre-WD lifetime is shorter, and usually 



Figure 10. Pre-WD lifetime in terms of the initial metallicity for 
white dwarf models with masses 0.50, 0.55, 0.60 and 0.65Mq. 

negligible, as compared to the cooling time on the white 
dwarf sequences once they evolved to lower luminosities. 

We can do a simple exercise by estimating the total 
age of a typical halo globular cluster. Kalirai et al. (2009) 
determined spectroscopically that the mass for the bright 
white dwarf stars in M4, at the beginning of the observed 
cooling sequence, is 0.592 ± 0.012Mq. Assuming that we 
do not know the metallicity of the cluster, we can estimate 
the age of the progenitor, and of the star, from the pre-WD 
ages listed in Table [T] The results are listed in Table (4] along 
with the estimated mass of the progenitor star. Note that all 
progenitor masses are low mass stars and do not experience 
many thermal pulses, between 1 and 5, during the AGB. 
The pre-WD age does not decrease monotonously, as we 
saw from Figs. [9land llOI For the metallicity value estimated 
spectroscopically for M4, Z ~ 0.001, the estimated values of 
progenitor mass for the bright white dwarfs is in agreement 
with the estimations of the stellar mass of the main sequence 
turn off, ~ 0.85M© (see Sec. 16.111 . 

6.1 The white dwarf cooling curve in globular 
clusters 

In Figure [TT] we show the total luminosity as a function of 
age during the white dwarf cooling sequences for four dif¬ 
ferent metallicity values Z — 0.0001,0.001,0.004 and 0.008. 
For each sequence we consider the total age of the star, i.e. 
the age measured from the ZAMS considering all previous 
stages to the white dwarf cooling curve, added to the cooling 
lifetime. The mass range for each metallicity is indicated in 
the plot. Finally, the horizontal dashed line represents the 
point in the cooling curve where crystallization in the central 
core begins, following of the Horowitz et al. (2010) prescrip- 
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Table 4. Age and progenitor mass for different metallicities for a 
white dwarf star with mass O.592±O.O12M0, estimated by Kalirai 
et al. (2009) for the bright white dwarf stars of M4. 


z 

Age [Gyr] 

MzaMs/Mq 

0.0001 

11.133 

0.83 

0.001 

10.587 

0.90 

0.004 

9.196 

0.97 

0.008 

8.834 

1.03 

0.01 

10.369 

1.04 

0.02 

7.499 

1.36 

0.04 

9.413 

1.22 


tion of the C—O phase diagram. The white dwarf sequences 
characterized by a progenitor metallicity of Z = 0.001 and 
Z = 0.01, originally computed in Renedo et al. (2010), were 
recalculated by considering the Horowitz et al. (2010) phase 
diagram. Note that with the Segretain & Chabrier (1993) 
phase diagram in the evolutionary computations, the effec¬ 
tive temperature at which crystallization begins at the cen¬ 
ter of the star is ~ 1000 K higher than that predicted by 
Horowitz et al. (2010) (see Romero et al. 2013 for details). 
Although white dwarf atmospheres at late stages of the cool¬ 
ing sequences are metal free, the final white dwarf evolution 
is not insensitive to the metallicity of the progenitor star. As 
we showed in previous sections, the pre-WD lifetime is larger 
for metal rich sequences as compared to low metallicity pro¬ 
genitors. In case of a IMq progenitor star, the progenitor 
lifetime goes from 5.9 Gyr to 11.1 Gyr when the initial metal¬ 
licity increases from Z = 0.0001 to Z = 0.01. On the other 
hand, as we showed in sec. EU the total amount of hydrogen 
left in the white dwarf remnant is usnally larger for lower 
metallicities, so a thicker hydrogen layer is expected for low 
Z sequences. 

Our lifetime computations during the white dwarf se¬ 
quences are consistent to those of BaSTI evolutionary code 
up to 2% at luminosities lower than log^L/Lo) ~ —1.5 
(Salaris et al. 2013). This difference is smaller than the un¬ 
certainties in cooling times attributable to the present uncer¬ 
tainties in the white dwarf chemical stratihcation. We also 
computed sequences with an initial mass IMq and metallic¬ 
ity of Z = 0.0001 and Z = 0.02 using MESA (Paxton et al. 
2011, 2013) and found lifetimes consistent with ours, with 
differences of ~ 0.1 Myr at 25 000 K on the cooling curve. 

HST observations of globular clusters in the Galactic 
Halo, M4 (Richer et al. 2004; Hansen et al. 2004), NGG 
6397 (Richer et al. 2006; Hansen et al. 2007) and 47 Tuc 
(Hansen at al. 2013), have uncovered several hundred mem¬ 
ber white dwarf stars. From a sample of galactic globular 
clusters, Krauss & Chaboyer (2003) derived an mean age 
of I 2 . 6 I 2 2 Gyr. More recently, Kalirai (2012) derived an 
age of the stellar halo near the position of the Sun to be 
11.4 ± 0.7 Gyr. Consequently, given these large ages, bright 
white dwarf stars are expected to have low stellar mass val¬ 
ues around 0.51 — O.SSMq (Renzini & Fusi Pecci 1988; Ren- 
zini et al. 1996). With the high precision data available we 
can employ the current observational constrains on age to 
derive some properties for the white dwarf sample in these 


systems. In particular, we used the total age determinations 
of three globular clusters to estimate the mass range on the 
white dwarf cooling curve and the minimum mass expected 
to have some fraction of its core in a crystallized state. 

The globular cluster NGG 6397 is one of the closest 
globular clusters (2.6 kpc) and is also a metal poor system 
with a iron to hydrogen ratio of [Fe/H]= —2.00 ± 0.01, cor¬ 
responding to Z = 0.0001. For this globular cluster Hansen 
et al. (2013) determined a total age of 11.7±0.3 Gyr. Con¬ 
sidering this total age we found that the white dwarf stars 
at high luminosities at the beginning of the cooling curve 
are characterized by a mean mass of 0.535 ± O.OOSMq cor¬ 
responding to an initial mass of ~ 0.85M@. On the other 
hand, the corresponding minimum stellar mass for a white 
dwarf to have undergone some degree of crystallization in 
the central core, as seen by Winget et al. (2009, 2010), is 
0.561 ± O.OO 7 M 0 , with a progenitor mass of ~ 0.95Mq. The 
observed white dwarf mass is not the same along the cool¬ 
ing curve. White dwarf stars with higher stellar mass cool 
faster, having lower luminosity values at a given age. In ad¬ 
dition, the pre-white dwarf lifetime decreases dramatically 
when the progenitor mass increases, so these stars reach the 
white dwarf stage in a shorter amount of time. 

For M4, Hansen et al. (2004) derived an age of 12.1 Gyr, 
with a 95% lower bound of 10.3 Gyr, by direct htting of the 
Color-Magnitude Diagram. Bedin et al. (2009) also deter¬ 
mine the cluster age of 11.6 ± 0.6 Gyr, consistent with the 
age from hts to the main-sequence turnoff (12.4±1.4 Gyr). 
The iron to hydrogen fraction for M4 is [Fe/H]= —l.lTO.Ol 
(Mucciarelli et al. 2011) corresponding to Z = 0.0013. We 
consider that the grid of sequences with Z=0.001 is repre¬ 
sentative of this cluster. By considering the age determina¬ 
tions from Bedin et al. (2009) we found that the mass at the 
top of cooling curve is around < 0.536 ± O.OISMq, in agree¬ 
ment with the mass determinations of Kalirai et al. (2009) of 
0.529 ± 0.012Mq . Note that the progenitor mass for a white 
dwarf with ~ O.536M0 is ~ O.92M0 at Z = 0.001. Also, 
those sequences with white dwarf mass ^ 0.553 ± O.OO 8 M 0 
have reached crystallization temperatures at an age ~ 11.6 
Gyr. So sequences with initial mass larger than IMq are 
already crystallized. The small difference between the age 
derived by Bedin et al. (2009) and that of Hansen et al. 
(2004) do not change considerably our results. 

The globular cluster 47 Tuc is more metal rich than 
both NGG 6397 and M4, with [Fe/H]=-0.75 (Carretta et al. 
2009), so we use a metallicity value of Z = 0.004 for this 
cluster. Hansen et al. (2013) determine a cluster age by ht¬ 
ting the properties of the cluster white dwarf population. 
They derived a total age of 9.9T0.7 Gyr at 95% conhdence. 
From our computations we found that the sequence charac¬ 
terized by a white dwarf mass of 0.537Mq has not reached 
the crystallization temperature (see Fig. [nj while all se¬ 
quences with white dwarf mass ^ 0.547Mq, and progenitor 
mass of ^ I.IMq, have a partially crystallized core at 10 
Gyr. Finally the white dwarf stars at high luminosity in the 
cooling curve are characterized by an stellar mass around 
0.524±0.005Mq corresponding to progenitors with O. 95 M 0 . 

Summarizing, the stellar mass values of the white dwarf 
stars at high luminosities derived for the three globular clus¬ 
ters are very similar, around ~ O. 53 M 0 , in agreement with 
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Figure 11. Total luminosity as a function of age for our set of sequences with four different metallicity of the progenitor star. The age 
corresponds to the total age of the model, by considering all the stages previous to the white dwarf cooling. The horizontal dashed line 
represents the point in the evolution where the crystallization process begins at the core, by means of the Horowitz et al. (2010) phase 
diagram for crystallization. The stellar mass values in solar mass units, from left to right, are: 0.519, 0.534, 0.550, 0.561, 0.569, 0.621, 
0.669, 0.708, 0.737 and 0.826 for Z = 0.0001, 0.505, 0.536, 0.553, 0.566, 0.593, 0.627 for Z = 0.001; 0.503, 0.515, 0.524, 0.533, 0.537, 
0.547, 0.553, 0.599 for Z = 0.004 and 0.498, 0.509, 0.516, 0.524, 0.531, 0.538, 0.552, 0.588 for Z = 0.008. 


previous estimates. Nevertheless, the progenitor mass at the 
ZAMS changes with metallicity. The progenitor mass for 
NGC 6397 is around 0.85Mq, while for M4 is ^ 0.92Mq. 
Since the total ages for these two systems are basically the 
same, within the uncertainties, the initial mass differences 
are a consequence of the different characteristic metallicities. 
47 Tuc is younger than NGC 6397 and M4 by ^ 1 — 2 Gyr 
and thus the mass of the main sequence is larger (O. 95 M 0 ), 
as expected. However the stellar mass at high luminosities in 
the white dwarf cooling sequences is the lowest of the sam¬ 
ple, indicating an enhanced mass loss rate, associated also 
to the higher metallicity. 


7 CONCLUSIONS 

In this work we studied the impact of metallicity on the 
theoretical white dwarf models, considering also its influ¬ 
ence on the previous stages of evolution through the analy¬ 
sis of the IFMR and the evolutionary lifetimes. To this end 
we compute a grid of full evolutionary sequences character¬ 
ized by metallicity values for the progenitor stars between 
Z — 0.0001 and 0.04, covering the Z range observed on the 
galactic disk and halo populations. We focused on stars with 
initial masses ranging from 0.85Mq to 3Mq, giving rise to 
white dwarf stellar masses in the range between ~ 0.5M© 
^ 0.87Mq. Our main results are the following: 
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• From our theoretical evolutionary computations, we 
found that the IFMR shows a dependence with the initial 
metallicity: the stellar mass of the remnant on the white 
dwarf colling curve decreases as the initial metallicity in¬ 
creases, indicating an strong influence of the metallicity on 
the mass loss rates during the RGB and the AGB phases. 
Our results are in good agreement with observations and 
semi empirical estimations of the IFMR, showing that the 
spread observed in the semi empirical data can be explain 
by means of the different metal abundances. 

Our analysis on the IFMR does not include the effect of 
extra-mixing episodes during the TP-AGB evolution. Then 
for most sequences dredge-up episodes, and its effects on the 
effective temperature, mass—loss rate, and core mass growth 
during the TP-AGB evolution, are not present. It follows 
that the dependence of the IFMR is explored mainly under 
the conditions of fixed initial metallicity, and photospheric 
G/0 ratio below unity during the TP-AGB phase, and thus 
no G-rich stars are formed from our models, except for se¬ 
quences with high initial mass and metallicity Z = 0.0001. 

• The mass of the helium core at the beginning of the 
central helium burning is smaller for more metal rich models. 
This will strongly influence the determination of the inner 
chemical composition of white dwarf stars with stellar mass 
~ 0.5Mq. Finally, the minimum mass for a star to undergo 
a helium flash decreases with metallicity. 

• We confirm the dependence with metallicity of the hy¬ 
drogen mass left at a white dwarf star, being in general 
thinner for models with higher Z progenitors. 

• We found that as the metallicity increases the pre-white 
dwarf lifetimes also increases, with age differences up to ~ 
5 — 6 Gyr for a ~ IMq star. Then, in order to estimate the 
age of a stellar population or a single object properly, we 
not only need to consider the total pre-white dwarf lifetime 
but also we must take into account its characteristic metal 
content. 

• Finally we employed our evolutionary computations to 
study three old globular clusters characterized with very 
different metallicity values: NGG 6397 {Z = 0.0001), M4 
{Z — 0.001) and 47 Tuc {Z = 0.004). By taking the age 
determinations from the literature, we estimate the mass of 
the white dwarf at high luminosities and the corresponding 
mass of the progenitor star, employing the IFMR with the 
corresponding Z value. We found that all three clusters have 
a similar white dwarf mass at high luminosities 0.53M©), 
in agreement with previous estimates. However, the initial 
mass values are different, being larger for 47 Tuc. This re¬ 
sult is an indication of the dependence of IFMR with the 
metallicity. 


In future works an extension to higher initial masses 
should be performed, from the theoretical point of view in 
order to study younger stellar populations, as open clusters. 
Also, this could be used to constrain the free parameters 
governing the theoretical mass loss rates that are currently 
used in stellar evolutionary computations. 
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